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ABSTRACT 

Context. Observations on 2011 August 9 of an X 6.9-class flare in active region (AR) 11263 by the Atmospheric Imaging Assembly 
(AIA) on-board the Solar Dynamics Observatory (SDO), were followed by a rare detection of vertical kink oscillations in a large- 
scale coronal active region plasma curtain in EUV coronal lines. The damped oscillations with periods in the range 8.8-14.9 min were 
detected and analyzed recently. 

Aims. Our aim is to study the generation and propagation of the MHD oscillations in the plasma curtain taking into account real¬ 
istic 3D magnetic and density structure of the curtain. We also aim at testing and improving coronal seismology for more accurate 
determination of the magnetic field than with standard method. 

Methods. We use the observed morphological and dynamical conditions, as well as plasma properties of the coronal curtain based 
on Differential Emission Measure (DEM) analysis to initialize a 3D MHD model of its vertical and transverse oscillations by imple¬ 
menting the impulsively excited velocity pulse mimicking the flare generated nonlinear fast magnetosonic propagating disturbance 
interacting with the curtain obliquely. The model is simplified by utilizing initial dipole magnetic field, isothermal energy equation, 
and gravitationally stratified density guided by observational parameters. 

Results. Using the 3D MHD model, we are able to reproduce the details of the vertical oscillations and study the process of their 
excitation by nonlinear fast magnetosonic pulse, propagation, and damping, finding agreement with the observations. 

Conclusions. We estimate the accuracy of simplified slab-based coronal seismology by comparing the determined magnetic field 
strength to actual values from the 3D MHD modeling results and demonstrate the importance of taking into account more realistic 
magnetic geometry and density for improving coronal seismology. 
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1. Introduction 

With the launch of the Solar Dynamics Observatory (SDO) the 
detection of waves in Extreme Ultraviolet (EUV) emission in the 
corona using Atmospheric Imaging Assembly (AIA) instrument 
became possible in many observed events with potential appli¬ 
cations to coronal seismology (CS) to understand the physical 
conditions in the corona. Using SDO/AIA EUV channels trans¬ 
verse coronal loop oscillations were observed to be ubiquitous 
in active region loops (McIntosh et al. 2011) and were analyzed 
in great detail with application to coronal seismology (e.g., As- 
chwanden & Schrijver 2011; Yuan & Nakariakov 2012; White & 
Verwichte 2012; Wang et al. 2012; Nistico et al. 2013; Verwichte 
et al. 2013b,a; Nistico et al. 2014). Nevertheless, vertical oscilla¬ 
tions in coronal loops or magnetic arcades are rarely observed in 
the EUV emissions in-spite of high-resolution and high-cadence 
recent observations by AIA, and past high-resolution observa¬ 
tions by the Transition Region and Coronal Explorer (TRACE) 
satellite. The vertical kink oscillations of a curved coronal loop 
are fundamentally different from the horizontal kink oscillations, 
since they are confined to the loop curvature plane, and can lead 
to the change of loop length. While horizontal kink oscillations 
were first detected by Aschwanden et al. (1999) and Nakariakov 
et al. (1999), and were studied extensively since then in many 
papers motivated by their potential use in MHD CS for local¬ 


ized coronal structures (e.g., reviews by Nakariakov & Verwichte 
2005; Andries et al. 2009; De Moortel & Nakariakov 2012; 
Liu & Ofman 2014, and references therein), only later Wang & 
Solanki (2004) have reported the first confirmed detection of ver¬ 
tical oscillations in a coronal loop using TRACE 195 A observa¬ 
tions of a coronal loop on the limb on 17 April 2002, with addi¬ 
tional cases analyzed in Wang et al. (2008). Vertical oscillations 
in a hot coronal loop following a flare/CME was reported for the 
first time by White et al. (2012) using SDO/AIA 131 A (~10 
MK) and 94 A (~6.3 MK) EUV bandpasses. Vertical oscilla¬ 
tions in cool coronal loops were also detected with SDO/AIA in 
conjunction with STEREO/EUVI. These oscillating loops were 
triangulated using STEREO/EUVI data, and oscillation proper¬ 
ties derived from SDO/AIA and STEREO/EUVI were utilized 
for refined coronal seismology (Aschwanden & Schrijver 2011). 
Recently, Srivastava & Goossens (2013) (hereafter, SGI3) have 
reported vertical oscillations in a coronal active region plasma 
curtain observed by SDO/AIA in cool coronal line of 171 A (~1 
MK). Vertical transverse oscillations were also observed recently 
in a coronal magnetic flux rope using SDO/AIA (Kim et al. 
2014). 

While several 2D numerical and analytical models of ver¬ 
tical kink oscillations in the coronal structures were published 
(e.g. Selwa et al. 2005, 2006; Verwichte et al. 2006; Gruszecki 
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et al. 2006; Selwa et al. 2007), the full three dimensional MHD 
modeling is required to study the complete interaction between 
the excitation mechanisms and the 3D mode polarizations of the 
oscillations in curved magnetic field geometry and the impul¬ 
sively excited propagating nonlinear fast magnetosonic waves 
launched by a flare that induces the oscillations. The 3D MHD 
modeling of a bipolar active region (AR) loop oscillations pro¬ 
vides further insights on the physical scenario and inherent cause 
for the somewhat rare detection of this phenomena (Selwa et al. 
2011b). At present, the application of coronal seismology is of¬ 
ten based on linear dispersion relation for waves in cylindrical 
slabs or otherwise simplified structures that compromises the ac¬ 
curacy of the method (e.g. Ofman 2009; Pascoe & De Moortel 
2014). In order to improve CS one needs to apply more realis¬ 
tic modeling of the oscillating structures. More realistic, bipolar 
magnetic field (e.g. Miyagoshi et al. 2004; McLaughlin & Of¬ 
man 2008; Selwa & Ofman 2009; Selwa et al. 2011a,b, 2013; 
Pascoe & De Moortel 2014) or extrapolated magnetogram based 
magnetic field (e.g., Ofman 2007; Schmidt & Ofman 2010) and 
gravitationally stratified density structures were modeled with 
3D MHD in order to study waves, clearly demonstrating the ad¬ 
vantage of considering 3D realistic magnetic field topology on 
unambiguous determination of seismologically derived coronal 
magnetic field. 

In the present article, we utilize recent observations by SGI3 
of vertical oscillations of a coronal curtain to setup the parame¬ 
ters of the 3D MHD model and study for the first time the exci¬ 
tation, propagation, and damping of the vertical oscillations in a 
coronal bipolar curtain with applications and quantitative tests of 
coronal seismology by comparing to 3D MHD modeling results. 
The paper is organized as follows: in Section 2, we present the 
brief summary of the observations, the numerical model is de¬ 
scribed in Section 3, the numerical results are in Section 4, and 
Section 5 is devoted to the Discussion and Conclusions. 

2. Summary of the Observations and DEM analysis 

Recently, the detection of a diffused and laminar plasma curtain 
at western equatorial corona that oscillates vertically during pas¬ 
sage of the coronal disturbances (EUV waves) triggered by X 6.9 
class solar flare (start: 07:48 UT; peak: 08:05 UT; end: 08:08 UT) 
in AR 11263 (N17 W69) in its vicinity on 9 August 2011 was 
reported and analyzed by SGI3 (cf.. Fig. la). The vertical kink 
oscillations are triggered near the flare energy release site within 
the curtain (8.8 min); in the middle both at the apex (14.9 min) 
and inside (13.3 min), as well as at the surface in the south-most 
part of the curtain (12.7 min). The properties of these excited 
waves strongly depend on the local plasma and magnetic field 
conditions of the plasma curtain. The right-panels of Fig. 1 show 
the temporal evolution of the velocity in the middle at the surface 
(Fig. lb: solid line), inner part (Fig. lb: dashed line), as well as 
at the surface in the southward part (Fig. lb, dotted-dashed line) 
of the curtain, and near the flare blast site (Fig. Ic). These veloc¬ 
ity oscillations in the plane of the sky are completely localized 
in the space at various parts of the curtain and derived here from 
the time-distance measurement data of SGI3. The animation of 
the oscillation observed with SDO/AIA 171A channel are en¬ 
closed online with the article. We carried out emission measure 
and temperature distribution analysis of the western equatorial 
corona and plasma curtain during 8-9 UT on 9 August 2011 us¬ 
ing the automated method developed by Aschwanden & Boerner 
(2011). Estimation of these plasma properties provided clues on 
the localized thermal conditions of plasma curtain. We use the 
full-disk SDO/AIA in all EUV channels at three times covering 


the whole oscillation period of the curtain, i.e., 08:07 UT, 08:28 
UT, 08:44 UT on 09 August 2011. We calibrate and clean the 
data using the aia_prep.pro subroutine of SSWIDL and co¬ 
align the AIA images as observed in various AIA filters by using 
the co-alignment test as described by Aschwanden & Boerner 
(2011). Six AIA EUV filter full-disk and co-aligned images (304 
A, 171 A, 193 A, 94 A, 335 A, and 211 A) covering tempera¬ 
ture in the 0.5-9.0 MK range, have been utilized to calculate the 
DEM and produce the temperature maps. The;^^ lies mostly be¬ 
low 1.0 in the most of the parts of curtain except its top laminar 
layer and background corona above the curtain. However, in the 
given limit of the uncertainties, within the curtain we have good 
DEM forward fit, and with the least error of estimation the mea¬ 
sure of temperatures and emission measures in the region of the 
EUV curtain. 

Fig. 2 shows that the plasma curtain (and its oscillatory parts) 
are nearly at isothermal temperature at the beginning of the oscil¬ 
lations (8:00 UT), and without significant temperature variability 
throughout the oscillation sequence. The temperature inside the 
curtain varies by about 15% - close to the observed uncertainty. 
It is evident that hotter material is located above the curtain - 
and this may explain the darkening of the 171A emission in this 
region seen in Fig. 1, since the peak response of this channel is 
~1MK (Lemen et al. 2012). The hot material in the line-of-sight 
may also lead to apparently increased DEM temperature inside 
the curtain. Inspecting the Emission Measure (EM) map propor¬ 
tional to the column density in the line of sight in Fig. 2 (left 
panel) shows that the height-dependence of the density is con¬ 
sistent with exponentially decreasing density with height. This 
observation justifies our isothermal and initially hydrostatic ap¬ 
proximation model of the plasma curtain that is valid to under¬ 
stand the evolution of the wave phenomena in the present study. 

3. The 3D MHD Numerical Model 

For the present study we adapt the numerical model of a bipo¬ 
lar AR initially developed by Ofman & Thompson (2002), and 
used extensively to study magnetosonic waves in AR loops (e.g. 
McLaughlin & Ofman 2008; Selwa & Ofman 2009; Selwa et al. 
2011a,b; Ofman et al. 2012a; Wang et al. 2013). Note, it is not 
the goal of the present study to model the detailed magnetic, 
density, or thermal structure of the curtain. By using this model 
we aim to improve the MHD seismology of the observed verti¬ 
cal kink oscillations using more realistic field and density struc¬ 
ture than used in CS based on wave dispersion from simplified 
slab geometry. The improvement is achieved by matching the 
parameters of the global bipolar structure of the curtains’ field 
and the gravitational density stratification to the observationally 
deduced parameters. The parameters of the model were adjusted 
such that the idealized model is in good agreement with the dy¬ 
namics of oscillating laminar and diffused large-scale plasma 
curtain due to the X 6.9 class flare generated disturbances re¬ 
ported by SGI3. In particular, the aspect ratio of the dipole was 
selected to fit the observed active region curtain aspect ratio, the 
temperature that determines the gravitational scale height of the 
density was selected to match the observed line emission tem¬ 
perature, the normalizing Alfven speed and the parameters of 
the velocity Gaussian pulse launched at the boundary were op¬ 
timized through parametric study to match the observed wave 
amplitudes - further details given below. 

In order to model the event, we solve numerically the 
time-dependent, nonlinear, finite-yS, resistive 3D MHD equations 
with gravity, using a Cartesian grid in the range (-3.5,3.5) x 
(-3.5, 3.5) X (1.0,6.0), in normalized units (the normalizations 
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Fig. 1. (a) The snapshot on 9 August 2011 at 08:07 displays the X 6.9 flare (red arrow), EUV wave (yellow arrow), plasma curtain (blue arrow), 
and disturbances passing in a narrow channel above it (green arrow) (adapted from SGI3). The temporal evolution of is available as a movie in the 
online edition, (b) The temporal evolution of the velocity in various parts of the plasma curtain as observed by SGI3. (c) The temporal evolution 
of the velocity near the flare site. 


EM Map 08:44 UT 



850 900 950 10001050 
X (orcsecs) 


23.0 

22.5 

22.0 



850 900 950 10001050 
X (arcsecs) 


Temperature Map 08:28 UT 



850 900 950 10001050 


X (arcsecs) 


Temperature Map 08:44 UT 



850 900 950 10001050 
X (arcsecs) 



Temperature Map 08:07 UT 


200 


Fig. 2. Temperatures maps obtained from the DEM analysis of the coronal curtain shown in Eig. 1 using six AIA EUV lines at times 08:07 UT, 
08:28 UT, 08:44 UT on 09 August 2011 (right three panels) The emission measure map at 08:44 UT on 09 August 2011 is shown in the left panel. 


are given below). The adopted method for integrating the partial 
differential equations is a modified Lax-Wendroff method with 
fourth-order stabilization term (see, Ofman & Thompson 2002). 
Furthermore, we use Powell (1994) method to correct for any 
numerical nonzero divergence of the magnetic field. The set of 
equations includes the contribution from gravity and the energy 
equation is eliminated by the isothermal approximation (poly¬ 
tropic index 7=1). While using 7 = 5/3 would be appropriate 
for thermodynamic equilibrium conditions not applicable in the 
corona, or when the coronal heating terms and losses are known 
in detail and computed in the energy equation, the polytropic 


(isothermal) approximation is based on the reasonable approxi¬ 
mation that the heating and losses terms are balanced exactly in 
the energy equation, and combined with the high coronal heat 
conduction lead to constant temperature. The nearly isothermal 
structure of coronal loops in quiescent active regions has been 
deduced from past observations (e.g., Aschwanden 2005). 

Moreover, observational evidence of coronal loop oscilla¬ 
tions shows that the effective polytropic index is close to unity 
(Van Doorsselaere et al. 2011), in line with our choice of poly¬ 
tropic index. This approximation is supported by the DEM anal¬ 
ysis of the present structure showing about 15% variability of 
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the temperature in the active region curtain in Fig. 2, on the or¬ 
der of the observational uncertainty. Since in the oscillation in 
the curtain are observed for about 40 min, the effects of resistive 
dissipation can be neglected. At the end of the oscillations the 
intensity structure of the curtain reverts mostly to the original 
form, indicating negligible net effect of the oscillations on the 
curtain. Non-isothermal (cooling) loop oscillations were studies 
in the past (e.g., Morton & Erdelyi 2009; Erdelyi et al. 2011; 
Ruderman 2011; Al-Ghafri et al. 2014). However, in the present 
study we neglect the effect of cooling, since the observed active 
region temperature remains nearly constants on the oscillation 
timescale, as evident in Fig. 2. 

Multiple processors (up to 144) are used to execute the code, 
which is parallelized using Message Passing Interface (MPI). We 
assume that the magnetic field perturbed by the flare can be rea¬ 
sonably described by perturbing the potential dipole (Fig. 3a) 
in the 3D computational domain, modeled on a 512x512x434 
grid. The high resolution resolves well the background density 
and magnetic field gradients as well as the propagating fast mag- 
netosonic pulse launched to model the effects of the flare, and 
minimizes the possible effect of numerical dissipation - this was 
tested by comparing to lower resolution runs. The foot-points are 
located at (+2.5, 0.0,1.0) in the box of normalized units, such ge¬ 
ometry allows to preserve the aspect ratio of the observed plasma 
curtain (~ 530 Mm length and height of ~ 100 Mm). 

The MHD equations are normalized using the following pa¬ 
rameters: length scale of a = R^/10, magnetic field at the foot- 
points Bq = 23 G, temperature Tq = 1 MK is the isothermal 
temperature in the model active region curtain, particle density 
no = 10^ cm~^ at Zmin = 1- These values are guided by the 
observationally determined values obtained by SGI3 and de¬ 
termine the Alfven speed of Vao = 15^6.5 km/s, the Alfven 
time ta = 44.1 5*, and the isothermal sound speed of Cg = 
(Po/Po)^^^ = 129 km/s. 

The observed oscillations shown in Fig. 1 and modeled here 
are for plasma emitting from Fe IX at 171A in EUV. Since the 
peak response of the 171A EUV filter is around IMK (Femen 
et al. 2012) the observed oscillating magnetic structure shows 
IMK plasma component in the multi thermal coronal plasma, and 
justifying the choice of temperature for the numerical model of 
the present observation. Further, as discussed above it is evident 
from the DEM analysis in Fig. 2 the DEM temperature is nearly 
constant inside the active region curtain during the observed os¬ 
cillation period. The selected value for the Alfven speed exceeds 
the value estimated by SGI3 using simplified (slab) coronal seis¬ 
mology by about 20%. We found, in fact, that we need to adopt a 
larger value by increasing the background magnetic field Bq, in 
order to match the observed oscillation amplitude. One must not 
forget that in our work we are modeling MHD oscillations of a 
complex curved magnetic structure such as a coronal curtain by 
using an idealized model for the AR based on a dipole magnetic 
field, thus it is expected that the model values of the physical pa¬ 
rameters improve on a simplified slab model, but may still differ 
from the actual values. 

The initial density profile is assumed to be gravitationally 
stratified and characterized by a normalized scale height H 

p=p,e^mo.z-z,,„)-um/H^ ( 1 ) 

where H = IkBToRg/ (lOGM^m^), Rg is the solar radius, ks is 
Boltzmann constant. To is the temperature, G is the universal 
gravitational constant. Mg is the solar mass and mp is the pro¬ 
ton mass. The initial density distribution is showed in Fig. 3a, 
as a cut in the vz plane (y = 0), where v is the direction of the 


initial injection of the idealized disturbance that models the ef¬ 
fects of the flare (see below). Some closed magnetic field lines 
that belong to the curtain in this plane are shown. Although the 
image in Fig. la shows higher emission in the curtain than in 
the corona above this is likely due to the temperature increase of 
the ambient plasma above the curtain as evident from the DEM 
analysis in Figure 2. The narrow band 171A AIA channel is sen¬ 
sitive to temperatures around IMK (Femen et al. 2012) and the 
hotter material evident above the curtain appears darker in this 
emission. However, inspecting the emission measure map of this 
region in Figure 2 (left panel) does not show sharp transition 
boundary and is consistent with exponentially decreasing den¬ 
sity used in our model. In our model the geometry of the curtain 
is modeled by the dipole field and the (nearly) frozen-in plasma 
carried by the fieldlines during the oscillations. Moreover, our 
analysis is focused on the regions inside the curtain with obser¬ 
vationally guided height obtained from the magnetic loops as¬ 
pect ratio and the possible small effect on the transverse waves 
in the interface region between the cooler and hotter region is 
not considered. 

In Fig. 3b, we show the normalized fast magnetosonic speed 

(in terms of Vao), V/ (x, z) = [Va (■^, z)^ + Cjj for perpendicu¬ 
lar propagation in the vz plane for y = 0, of the initial state. As it 
is clear from the image, Vf speed decreases extremely fast with 
the distance from the footpoints, reducing to 25% of the maxi¬ 
mum value at about just half the height of the grid (z = 3.5). This 
is because, for Cg Va, as we can reasonably assume in this 
case, Vf behaves basically as Va, thus decreasing SiS B ^ 1/r^ (if 
far enough from the dipole) and being inversely proportional to 
which, in turn, decreases nearly exponentially with z and is 
characterized by a scale height H that depends upon the obser¬ 
vationally guided choice of the initial temperature of the coronal 
plasma in the curtain. The location of the points LI - L4 used 
for temporal evolution diagnostic are marked with ’+’ symbols 
in the xz plane, where the 3D coordinates in normalized units of 
the points are LI = (-1.75,-3.35,1.15), L2 = (0.0,1.7,1.3), 
L3 = (0.0,1.7,1.85X L4 = (1.75,-0.9,1.65). The choice of 
these points is consistent with the observing points in SGI3 
shown in Fig. 1. 

Boundary conditions are imposed on each surface at the 
edges of the 3D computational box. Open boundaries with zero- 
order extrapolation for every plasma variable are placed at all 
boundaries except for the lower boundary (z = Zmm), where 
we impose a condition of wave reflection simulating line-tied 
boundary, and for the interface between the integration domain 
and the external disturbances. Since modeling the flare is beyond 
the scope of the present study in order to trigger the vertical kink 
oscillations in the plasma curtain as observed by SGI3, we per¬ 
turb the AR initial state with a Gaussian shaped velocity pulse 
applied as time dependent boundary condition for short duration 
compared to typical normal mode periods (similar in form to the 
ones used in Selwa & Ofman (2009); Selwa et al. (2011a)) in the 
V direction (see Figure 3a) 


Vx = Ayexp 


i Z Zmin 

\ 


exp 


?-(?! +2(5r)V 

Jt / 


( 2 ) 


applied at the yz plane at v = Xmin for the duration of the pulse 
St = t2 - ti = It A- The pulse magnitude is normalized in units 
of Vao, and the parameters are Ay = 0.18, = 0.7, where L is 

the pulse start time at 0.5ta after the model initiation time t = f). 
The value of provides the width of the Gaussian pulse in z and 
deposits the energy in the lower part of the active region curtain, 
consistent with observations. The values of the pulse parameters 
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Fig. 3 , (a) Idealized model of the AR as a magnetic dipole. The red box shows the boundary plane at x = 0 impacted by the propagating 
disturbance that models the effect of the flare on the active region curtain. The red arrow shows the direction of the impacting velocity pulse applied 
as the time-dependent boundary condition on this plane. The orientation of the model is turned by 90^ CCW with respect to the observation shown 
in Fig. 1. Bottom: Normalized initial density distribution overlayed with some of the closed field lines that belong to the curtain are shown, (b) and 
initial fast magnetosonic speed (c)m the xz plane at y=0. The blue (white) ‘-1-’ symbols mark the locations of the points LI - L4 used for temporal 
evolution diagnostic (see below). 


have been chosen by performing a parametric study that best fits 
the observed wave amplitude at point L1-L2. For example, we 
have repeated the same model runs with several valued of Ay in 
the range Ay = 0.05 to Ay = 0.25, and found through interval 
halving of this parameter that Ay = 0.18 provides oscillation 
amplitudes that agree well with observations. However, the pe¬ 
riodicity and the damping rate of the waves are not sensitive to 
the choice of the pulse parameters for small 6t compared to the 
periods. The value of were varied in several runs in the range 
19-26 G, and as expected, it was found that the oscillation period 
is proportional to Bq with other parameters fixed, providing best 
fit for ^0 = 23 G. 

4. Numerical Results 

In this section we highlight the results of our MHD model of 
flare-induced vertical kink oscillations excited by means of a lat¬ 
eral Gaussian-shaped pulse, modeling the effects of the kinetic 
energy release by the X 6.9-class solar flare. 

Fig. 4 shows the evolved perturbation in terms of plasma 
velocity, perturbed plasma density (defined as Ap = [p{t) - 
p(0)]/p(0), and the magnetic field in the xz plane. The snapshot 
has been captured at r = 14.7 ta, at y = 0, where x is the in¬ 
jection direction of the velocity pulse. The form of the velocity 


pulse and the associated compression of the density are shown 
in Fig. 4a-b. The current density associated with the fast magne- 
tosinic wave front produced by the pulse, and the corresponding 
changes of the magnetic field direction are evident in Fig. 4c-d, 
respectively. The response of the curtain produces velocity in the 
negative x direction due to reflection and the restoring Lorentz 
force of the perturbed dipole field (see. Fig. 4a) with smaller 
magnitude than the pulse. The Gaussian pulse launched at the 
boundary causes rarefaction of the adjacent density region and 
propagation of nonlinear (steepened) compressive fast magne¬ 
tosonic wave that perturbes nonlinearly the magnetic structure 
of the curtain inducing oscillations. The propagation and the re¬ 
flection of the wave is nonuniform due to the bipolar structure 
of the magnetic field and the corresponding structure of the fast 
magnetosonic speed (shown in Figure 3b). As the fast magne¬ 
tosonic wave leaves the computational domain the oscillations 
of the curtain proceed to relax approaching normal mode, and 
eventually damp as the transverse waves leak out of the struc¬ 
ture, and the final density and magnetic field structure returns 
close to the initial structure. 

In Fig. 5, we display the temporal evolution of vertical fluctu¬ 
ations in the vertical plasma velocity component at locations 
LI - L4, very similar to what has been shown by observations 
(Fig. 1). We observe variations in the wave phase speed at differ- 
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Fig. 4. Snapshot of the cuts in the xz plane at the center of the curtain obtained with the 3D MHD model Bit - 14.7ta showing the propagating fast 
magnetosonic disturbance produced by the boundary pulse that models the effects of the flare, (a) The perturbed density Aplp(t = 0). Fieldlines 
associated with the active region curtain are shown. The distortion of the held lines by the propagating fast magnetosonic disturbance is evident, 
(b) The velocity magnitude and the direction (arrows are shown for |v| greater than 6% of the absolute maximal velocity), (c) The current density 
map f. (d) The magnitude of the magnetic held and the direction. The results show the effects of the initial induced wave pulse and the response 
of the curtain. 


ent points, not necessarily aligned with the y axis. The reason of 
this choice of ‘observing points’ in the model comes from direct 
observations of the event by SGI 3, suggesting that the bundle of 
loops experiencing the vertical kink oscillations might be quite 
tilted to the bottom boundary surface, and the oscillations are not 
necessarily co-planar. 

The dash-dotted curve in Fig. 5 shows vertical transverse 
kink oscillations near the boundary plane where the flare pro¬ 
duced pulse is modeled, i.e., nearest to the flaring site (LI) and 
taking place in the deeper layers of the plasma curtain (z = 1.3, 
~ 20 Mm above the photosphere). This oscillation has been 
interpreted by SGI3 as the first overtone with a period of 8.9 
min (cf.. Fig. Ic for velocity time dependence). For oscillations 
in the northward side of the plasma curtain, our MHD model 
yields a fundamental period Pq = In/co = 9.2 min with 99.99% 
confidence level, in surprisingly good agreement with SDO/AIA 
observations, given the simplicity of the present magnetic field 
model. 

The solid curve in Fig. 5 shows the plasma vertical velocity 
versus time at locations L2 and the dashed curved shows L3: near 


the loop apex and mid-way from the flaring site (y = 0), respec¬ 
tively. In this region there is observational evidence for oscilla¬ 
tions in both deep and surface layers of the coronal curtain (~ 20 
Mm and ~ 60 Mm over the surface). According to SGI3, the 
large-scale disturbance at depth could represent the fundamen¬ 
tal kink mode (but, possibly dominated by the first harmonic, 
see below) with a period of 13.3 min, while the non-decaying 
surface disturbance perturbs only the magnetic upper sheet of 
the curtain and is characterized by a period of 14.9 min. Once 
again our MHD model of the event seems to reproduce the am¬ 
plitude and period of the perturbations with a very good agree¬ 
ment: A=34.8 km/s, Pq = 13.5 min for the deep mode; A=40.1 
km/s, Pq = 14.8 min for the surface mode, where the periods are 
obtained from the periodorgram analysis (Scargle 1982), and the 
amplitude is obtained from sinusoidal least-square-flt with close 
valued of the periods. The time-profile of associated observed 
velocity perturbations can be seen in Fig. lb. 

Lastly, we show perturbations in the top layers of the coro¬ 
nal curtain, L4, shown with short dashes in Fig. 5, southward and 
opposite to the flaring site (y = 1.75). Vertical kink modes have 
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Fig. 5. Temporal evolution of the model velocity fluctuations in the yz 
plane at different locations along the x axis. LI: near flare blast deep 
layer oscillations (-1.75, -3.35,1.15); L2: near apex deep layer oscil¬ 
lations (0.0,1.7,1.3); L3: near apex surface oscillations (0.0,1.7,1.85); 
L4: southward surface oscillations (1.75, -0.9,1.65). 



Fig. 6. Temporal evolution of the model velocity oscillation in the yz 
plane at L2 (’-h’ symbols) fitted with a single frequency sinusoidal func¬ 
tion (green) and the damped sinusoidal function (blue). The resulting 
period for the damped sinusoidal is 12.8+0.8 min and the best-fit damp¬ 
ing time is in the range 6.2-11.7 min. The best-fit period is 12.4 + 0.8 
min for the sinusoidal fit. 

been observed in this region as well, with a period of 12.7 min 
by SGI3. Analysis of our modeled oscillations at this location 
yields A=46.1 km s“^ and Pq = 11.7 min, in good agreement 
with observational measurements (cf., Fig. lb). The detection of 
the fundamental and second harmonic in the oscillations is im¬ 
portant since their ratio in a coronal loop can be used for the de¬ 
termination of the gravitational scale height of the density (e.g., 
Andries et al. 2005; Ruderman et al. 2008; Verth et al. 2008). 
This effect is expected to be significant in long stratified coronal 
loop structures. However, in the present study we deal with os¬ 
cillations of a bipolar structure in short magnetic field line loops 
(see, Fig. 3), therefore, the magnitude of the second harmonic is 
found to be weak. Guided by SGI3 observations, high-contrast 
density structure of individual loops is not included in the model 
of the curtain. 

In Figure 6 we show the results of least-square exponentially 
damped single-frequency sinusoidal fit to the modeled oscilla¬ 



Period [min] 

Fig. 7. The results of the periodogram analysis (Scargle 1982) of the 
oscillations at L2 (solid) and L3 (dashes). The confidence level of the 
dominant period is shown with the dash line for each periodogram. 

tion at L2. The fitted function is of the form 

HO = sin{(jjt + 0), (3) 

where A is the damping time and oj = InjP is the oscillation fre¬ 
quency. We also show an undamped single-frequency sinusoidal 
fit. The damped sinusoidal fits produces a period of 12.8+0.8 
min and the damping time in the range 6.2-11.7 min. The fit pro¬ 
duces an amplitude of 355 km s“^ and the phase 0 = 4.6 Rad 
at r = 0. However, the values of the amplitude and the phase 
at r = 0 are not directly relevant in this analysis, since the os¬ 
cillations at L2 start at r ~ 9 min due to the travel time of the 
pulse from the excitation boundary to L2. Due to the nonlin¬ 
earity in the initial stages of the oscillations that produces large 
initial velocity fluctuation, the fitted exponential damping time is 
overestimated. However, it is interesting to note that the damp¬ 
ing time is of the same order as the period of the oscillations. 
Since the large amplitude of the initial nonlinear pulse may af¬ 
fect the damped sinusoidal fit, we have also fitted undamped si¬ 
nusoidal to the modeled velocity oscillations. At point L2 we get 
a period of 12.4+0.8 min for the undamped sinusoidal. We note 
that the values of the periods obtained from the damped sinu¬ 
soidal, undamped sinusoidal, and the periodogram analysis are 
close within the respective error bars. 

The periodogram analysis (Scargle 1982) of the modeled os¬ 
cillation time sequences for points L2 and L3 are shown in Fig. 7. 
It is evident that the oscillations are dominated by a single fre¬ 
quency with similar results at points LI and L4 (not shown). 
While there is some evidence of second and higher harmonics 
in the analysis, the statistical confidence level of these harmon¬ 
ics is low < Icr. The low statistical significance of the second 
harmonic may explain why the second harmonic is rarely de¬ 
tected in observations of loop oscillations. We have also per¬ 
formed curve fitting of single-period sinusoidal and exponen¬ 
tially damped sinusoidal function to the oscillations and found 
that the values of the fitted periods are in excellent agreement 
with the periodogram values. We And that the damping time of 
the oscillations is on the order of the oscillation period. 

Table 1 summarizes our results from the model velocity os¬ 
cillation, and compares them to SGI3 results obtained using 
SDO/AIA data. The modeled error bar of the periods obtained 
from least-square fitting is below 3%. However, since there is no 
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statistical variability in a computational model, the meaning of 
error bars is ill-defined in this context. Therefore, the computa¬ 
tional values are shown without error bars. The last two columns 
of Table 1 show the comparison between the mean magnetic field 
along the magnetic field lines (loops) that pass through the mea¬ 
surement points in the 3D MHD model (Bmhd), and the mag¬ 
netic field inferred from application of coronal seismology (CS), 
Bcs • Since the gravitationally stratified density and the magnetic 
field vary significantly along the loops (e.g., Ofman et al. 2012a) 
we use the average values along the magnetic loops. We are 
comparing the values of the average magnetic field along the 
magnetic loops that pass through the points L1-L4 in our MHD 
model, with the average value of the magnetic field determined 
from coronal seismology that uses a slab model for these loops 
(Equation 4). Since the density along a loop is close to the sur¬ 
rounding density we use Ck ~ Va in the low-yS plasma, and we 
calculate the loop-averaged density to find Bcs • Since the wave¬ 
length of the fundamental mode oscillations is twice the loop 
length, the waves cannot be used to determine magnetic field 
value at a particular point of the variable magnetic field in the 
loop (see. Figure 3 in Ofman et al. (2012a), for an example of 
magnetic field, density, and fast magnetosonic speed variability 
along a loop in the bipolar field model). The long wavelength 
oscillations provide information only on the average magnetic 
field in such loops, and these value are given in Table 1. 

The event animation from SDO/AIA shows how the plasma 
curtain oscillations are combinations of the fundamental and sec¬ 
ond normal modes. However, the latter contribution seems to be 
predominant at deep and surface levels. Therefore, CS predic¬ 
tions are computed using the expression for second harmonic 
oscillations of the ideal vertical kink mode without density con¬ 
trast across the loop (consistent with observations that show little 
density contrast), and average values along the loops 

Bcs = T, ( 4 ) 

where (p) is the average density of the loop, L is the loop length, 
Pq is the oscillation period in the MHD model, and Bqs is the 
average magnetic field magnitude along the loop. 

The periods of the observed and modeled oscillations show 
good agreement. Comparing the MHD and coronal seismology 
predictions of the mean magnetic field, we conclude that the dif¬ 
ferences are within 39% - 47% at 3 out of 4 observation loca¬ 
tions (Fig. 5b-d and row 2-4 in Table 1). This level of agreement 
was also found by Aschwanden & Schrijver (2011) when com¬ 
paring Bcs values to extrapolated magnetogram field values and 
by Pascoe & De Moortel (2014) in 3D MHD arcade model test 
of CS. The poorer agreement between the values near the Flare 
Blast location (Fig. 5a and row 1 in Table 1, Bcs I Bmhd ~ 4) is 
likely due to two factors: the large tilt of the magnetic loop that 
passes through the observation point to the yz plane, causing the 
loop to interact with the lower boundary surface, and the possi¬ 
bility that the observed oscillation is not a normal harmonic but 
rather a driven propagating mode due to the close proximity of 
the boundary that launches the Gaussian velocity pulse. 

5. Discussion and Conciusions 

The article demonstrates that an ideal dipole approximation of 
the magnetic field geometry is a good approximation of the mag¬ 
netic field structure in an AR plasma curtain on a large spa¬ 
tial scale over the western equatorial corona observed on 2011 
August 9 by SDO/AIA. This large-scale magneto-plasma sys¬ 


tem, when encountered with the flare generated propagating dis¬ 
turbances, produces multiple harmonics of vertical kink oscil¬ 
lations in the different parts of the curtain, as well as the sur¬ 
face oscillations. The 3D MHD numerical simulation domain 
includes a magneto-plasma system with gravitationally stratified 
density, and its interaction with the velocity pulse introduced at 
the boundary mimicking the flare disturbance, which propagates 
into the modeled AR curtain and produces magnetic field dis¬ 
placement and velocity oscillations of vertical kink waves. The 
results of the numerical simulation match well with the obser¬ 
vations of SGI3, while it is initialized with parameters based on 
the observed conditions, and thus improve the understanding of 
the wave activity. 

While the present 3D MHD model produces overall evolu¬ 
tion of the vertical transverse waves in good agreement with ob¬ 
servations - the scope of the model is limited, and is not intended 
to study flares and the detailed heating and cooling dynamics 
of an active regions. The main limitation of the model is the 
use of dipole field that is only approximately applicable to this 
observation, and this could be improved by using extrapolated 
magnetogram field (e.g. Ofman 2007) in future studies. The de¬ 
tailed effects of coronal heating, thermal conduction, and radia¬ 
tive cooling are not modeled and only approximately accounted. 
The use of polytropic energy equation with y = 1 is only appli¬ 
cable to approximately isothermal conditions where the above 
effects nearly balance and the heat conduction is high along the 
field. The initially hydrostatic density structure while providing 
overall reasonable approximation, does not capture the full com¬ 
plexity of active region and individual loop density structure and 
neglects the potentially important role of upfiows and downflows 
that can affect the wave dynamics. 

The damping of the oscillations in the model is due to refrac¬ 
tion and leakage of the fast mode wave flux as demonstrated in 
previous work (Ofman & Thompson 2002; McLaughlin & Of¬ 
man 2008; Selwa & Ofman 2009; Selwa et al. 2011a) since non¬ 
ideal damping processes, such as resistive or numerical dissipa¬ 
tion are small in the present model. It is evident that the observed 
oscillations appear to damp on a timescale similar to the modeled 
evolution. The effects of mode conversion (i.e., resonant absorp¬ 
tion) as discussed by SGI3 are possible but not strongly sup¬ 
ported since the observed AR curtain shows little evidence of 
individual loop structures (with high density contrast) necessary 
for this damping mechanism. The modeled low-jS plasma curtain 
contains strong magnetic field and the density is gravitationally 
stratified to mimic the observations and to best fit with the real¬ 
istic solar atmosphere. Therefore, a significant fast-mode speed 
gradient is set across the plasma curtain from its deeper layers 
towards the apex. This may set in a possibility of internal re¬ 
fraction of the waves and transfer of the fast mode (kink) wave 
energy towards the surface and outside the curtain (as demon¬ 
strated numerically by, e.g., Ofman & Thompson 2002). Some¬ 
what alternative interpretation of the refraction process can be 
found in Goossens et al. (2002, 2011, 2012), adapted by SG13 
initial analysis. 

The 3D MHD model of the waves in the coronal curtain al¬ 
lows us to test the application of simplified coronal seismology 
as it is widely used (see the reviews, e.g., Nakariakov & Ver- 
wichte 2005; Liu & Ofman 2014) to the present model and com¬ 
pare the results of CS with the actual values of the magnetic 
field used in the 3D MHD model. We find that simplified coro¬ 
nal seismology provides magnetic field values within 40%-50% 
of the actual values. Noting, that the magnetic energy of AR im¬ 
portant for Space Weather impact analysis scales as B^, this dis¬ 
crepancy amount to factor of ~2 of the magnetic energy content. 
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Location of observation 

Observed periods [min] 

MHD period [min] 

Bmhd (G) 

Bcs (G) 

Flare Blast (LI) 

8.9 

9.2 

2.2 

8.4 

Apex Inside (L2) 

13.3 

13.5* 

3.1 

4.3 

Apex Surface (L3) 

14.9 

14.8 

2.9 

4.3 

Southward Surface (L4) 

12.7 

11.7 

3.4 

4.9 


Table 1. Comparison between observed and model oscillation periods for the four detected location in SGI3. The values of the MHD periods 
were obtained using periodogram analysis (Scargle 1982) and the period with the highest intensity was chosen as the proxy for the normal mode. 
The penultimate column shows the magnetic field strength obtained from the 3D MHD model, and the last column is from simplified coronal 
seismology (CS). *At point L2 the average of two periods near the peak were taken since they have close high confidence level (see, Fig. 7). 


This result is consistent with previous tests of CS in more simpli¬ 
fied setup (e.g., Pascoe & De Moortel 2014, in magnetic arcade 
without gravity) and with comparisons to extrapolated magne¬ 
togram field (Aschwanden & Schrijver 2011). Thus, in order to 
improve CS and provide more accurate values of the magnetic 
field in ARs it is necessary to consider the properties of waves 
in complex, curved, and more realistic magnetic structures with 
realistic global density structure as demonstrated in the present 
study. 
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